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The optimized effective potential (OEP) method allows for 
calculation of the local, effective single particle potential of 
density functional theory for explicitly orbital-dependent ap- 
proximations to the exchange-correlation energy functional. 
In the present work the OEP method is used together with 
the approximation due to Krieger, Li and Iafrate (KLI). We 
present the first application of this method to polymers. KLI 
calculations have been performed for the insulating polyethy- 
lene and the results have been compared to those from other 
orbital-dependent potentials. Various properties of the band 
structure are also calculated. The single-particle band gap 
strongly depends on the basis set with larger basis sets yield- 
ing narrow gaps. For certain physical quantities such as the 
total energy and the exchange energy, the various orbital- 
dependent exchange-only and Hartree-Fock results differ only 
slightly. For the highest occupied orbital energy the differ- 
ence is more pronounced. In order to get the right band gap 
in OEP the exchange contribution to the derivative discon- 
tinuity is calculated and added to the Kohn-Sham gap. The 
corrected gap obtained by the KLI approach is 12.8 eV com- 
pared with the Hartree-Fock and experimental values of 16.6 
and 8.8 eV, respectively. We observe, however, the strong 
dependence of the derivative discontinuity on the basis set. 



1. Introduction 

The optimized effective potential (OEP) method was 
invented by Talman and Shadwick [Q , following the orig- 
inal idea of Sharp and Horton In the now classic 
OEP approach a given total energy functional depend- 
ing explicitly on single-particle orbitals is minimized un- 
der the constraint that these orbitals are solutions of a 
single-particle Schrodinger equation with a local, effective 
potential. The exchange-only OEP can be interpreted 
as the implementation of the exact exchange-only den- 
sity funtional theory (DFT) ||§. Very recently OEP is 
extended to include electron correlation effects via per- 
turbation theory |^,||. In this theory both the correla- 
tion energy and the potential are expanded using stan- 
dard perturbation theory. In this way a self-consistent 
parameter-free KS ah initio scheme can be given. The 
exact energy functional and the exchange-correlation po- 
tential provided by the Hohenberg-Kohn theorem can be 
calculated order-by-order by means of the OEP pertur- 
bation theory Q. 

To the best of our knowledge, full OEP calculations 



have so far only been performed for atoms and for 
semiconductors [|J at the exchange-only level of the the- 
ory. The incorporation of electron correlation effects into 
the OEP formalism has been the subject of few studies 
[^,0. Most of these studies, however, have used the 
complicated self-energy formalism. 

To reduce the numerical complexity of the OEP 
method, Krieger, Li and Iafrate (KLI) introduced a re- 
markably accurate analytical approximation to OEP and 
transformed the original OEP equations into an easily 
manageable form . Applications to atoms and to solids 
show that the KLI method yields results that are nearly 
identical to those of OEP The KLI method has 

been successfully applied for molecules |l2| and for cer- 
tain semiconductors Jl||14| using pseudopotentials. 

In our previous study we have carried out calculations 
by means of various orbital-dependent exchange poten- 
tials which are accurate approximations of OEP and re- 
quire less computational effort fll5| . We found that the 
various orbital dependent potentials provide a narrow 
DFT gap for polyethylene (PE). KS theory will typically 
give too narrow valence bands and too compressed con- 
duction bands [|llj . The KS theory will only give an exact 
result for the highest occupied energy eigenvalue, there- 
fore, calculated states will in general not have the desired 
accuracy as one might hope In addition to this, it 

has been shown that the difference between the highest 
occupied and lowest unoccupied DFT eigenvalues in the 
iV-electron system (which is what is usually called the 
DFT band gap) is not the true quasiparticle band gap, 
but differs from it by the discontinuity in the exchange- 
correlation potential when an electron is added to the 
system 0,^3). Therefore, the so-called derivative discon- 
tinuity corrections must be taken into account in order to 
get the right gap p 15 2|]. Consequently, it is desirable 
to calculate the band structure using the KS method, and 
then simply add on a single-shot correction to all states 
]IT| , p^[ . In this way much better agreement is found with 
the Hartree-Fock gap when the exchange corrections are 
considered to the KS eigenvalue band gap |l5| . None of 
the state-of-the-art exchange-correlation functionals de- 
pending only explicitly on the density show the above- 
mentioned discontinuity in the exchange-correlation po- 
tential and therefore do not yield the correct values for 
the band gaps of the insulating solids jll[] . 

In the last decades quite a number of publications 
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appeared on polymer electronic structure using mainly 
Hartree-Fock ab initio methods [0 . Results are reported 
only in recent years, however, using DFT methods like 
local density approximation (LDA) [E6j. The general- 
ized gradient approximation (GGA) ]2l| resulted in only 
modest improvement over LDA gap in semiconductors 
& 

In this paper, we briefly review the theory of OEP and 
KLI methods. This work can be taken as the systematic 
continuation of our previous work on polymer electronic 
structure obtained by various orbital dependent local ex- 
change potentials |15|. Particular attention is paid to 



the correction of the OEP eigenvalue band gap. We will 
present the first results obtained by the KLI orbital de- 
pendent exchange potential for polymers. We give the 
calculated band gaps, band widths, band structures and 
total energy of polyethylene within the Kohn-Sham den- 
sity functional scheme using various orbital dependent 
exchange potentials. 



2. Basic Formalism 

The following spin-restricted Kohn-Sham type mean- 
field approach is employed for polymers (the extension of 
the theory for spinpolarized problems is straightforward) : 
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V 2 + v s (r)]ui(r) = ejUi(r) 



(1) 



where {ui(r)} are the single particle orbitals and v s (r) is 
the Kohn-Sham effective single particle potential fl . 



v s (r) = v ext (v) + v H (r) + v xc {r) 



(2) 



The Kohn-Sham one-body potential contains the exter- 
nal, Hartree and exchange-correlation potential. v ext (r) 
represents the Coulomb potential of the nuclei. vh(j) is 
the classical Coulomb repulsion of the electron charge 
density.The exchange-correlation potential v xc (r) is a 
nonclassical term. It is formally defined as functional 
derivative of the exchange-correlation energy E xc [p\: 



=(r,[d) = 



SE xc [p] 
Sp{r) 



(3) 



We seek an expression which contains the exact exchange 
potential and correlation can be added easily on top of 
the exact exchange-only theory using approximate cor- 
relation energy functionals J^]. It must be emphasized, 
however that the correlation energy functionals either in 
their local or gradient corrected form may not work in 
combination with exact exchange. This is because basi- 
cally all common correlation functionals do not have a 
long-range component in the corresponding correlation 
hole (the combined exchange-correlation hole is typically 
short-ranged). The Colle-Salvetti gradient corrected cor- 
relation functional p8[, which provided excellent results 



for atoms, performed rather badly for molecules for the 
abovementioned reason |^Jl2]|. Approximate correlation 
functionals which are derived from the homogeneous 
or inhomogeneous electron gas model or obtained non- 
empirically from sum-rule conditions Jl6[ | have an incor- 
rect long-range tail (only dynamical correlation is ac- 
counted for). The correlation energy functionals provide 
improper local behaviour in the bond midpoint regions 
in molecules Jl6|,[l9],^0| which is mainly due to the lack 
of non-dynamical correlation effects fl6|| . Therefore, in 
those methods which work with exact exchange (HF and 
OEP) the correlation effect can only be accounted for 
in a more complicated way than in the so-called second 
generation of DFT 

In the optimized effective potential method (OEP) in 
contrast to the ordinary DFT, the exchange-correlation 
energy is approximated as an explicit functional of or- 
bitals and only implicit functional of the electronic den- 
sity O]. The starting point of the OEP method is the 
total energy functional 

S t ° f P [p] = f:/rfr<(r)(-iv 2 ) Ul (r) 

dvp(v)v ext {T) + \ J dvdv'&^l+E^liui}]. (4) 

We make use of the chain rule for functional derivatives 
to obtain from Eq. (||) 



OCC r, 

^ P (r,W) = E/^ 

i—1 



,6 E° EP [{ Ui }] 5 Ui (r') 
5ui(r') 5p(r) 



+ c.c. (5) 



where formally are the orbitals which are, how- 

ever, implicit functionals of the density. Applying the 
functional chain rule again and after some algebra one 
obtains the following integral equation: 

OCC p 

£>i(r) / dr'K(r')-^(r')]G s ,(r,r'K(r') + c. C . = 

i—1 

(6) 

where 



and 



1 SEg EP [u t ] 
<(r) Ui (r) 

u fc (r)w*(r') 



£i - £fc 



(7) 



(8) 



The integral in Eq. (^|) is the fundamental expression 
for v xc in OEP and can be shown to be equivalent 
to Eq. (||) Q. There is no known analytic solution 
for «a; C [{ui}]. Therefore, only numerical solutions are 
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available for spherical atoms [p]J^,|7|, p9| , ^0t and for solids 
P, fU3|fil| ] using atomiclike wave functions and assuming 
integer atomic occupation numbers within each muffin- 
tin sphere (linear augmented plane- wave method) [pj. 
These numerical solutions of Eq. (|^) are confined to the 
exchange-only OEP and correlation has been taken into 
account only via approximate local functionals B. In 
the work of Krieger and co-workers the OEP integral 
equation is analyzed and a simple approximation is made 
which reduces the complexity of the original OEP equa- 
tion significantly and at the same time keeps many of 
the essential properties of OEP unchanged 0. Krieger, 
Li and Iafrate gave an exact expression transforming the 
OEP integral equation (j^) into a manageable form. They 
obtained the following, still exact expression for v xc (r): 



vZ EP (r)=v s xc (r)+ J2 



Mr) ( _ ; 



OEP — 



lgV[pi(r)V Ui (r)] 



Vi) + 



(9) 



p(r) 



dYp i (r)v l {r), 



where m — occ is the highest occupied one electron en- 
ergy level, pi — \ui\ 2 is the partial density, uf (r) is 
the Slater's potential (its correlation part is not known 
exactly) 



«f(r) 



1 



2p(r) 



A|2 



dr 



/l7(r,r')| 



(10) 



where 7(r, r ) is the first-order density matrix. The sum- 
mation runs over the orbital index for all the occupied 
orbitals up to the highest occupied mth orbital (Fermi 
level) . The function pi is defined by 
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dr>°*V)-«i(rO]G«(ryK(r'). 

(11) 



In practical applications the last term in Eq. (9) turned 
out to be quite small in atomic systems and has small 
effect only on the atomic shell boundaries [0. This ad- 
ditional term's average over the mth orbitals is zero . 
Neglecting this term one obtains the exchange-only KLI- 
approximation which, after some algebra, can be written 
in the following form |tJ 

occ — 1 I / occ— 1 

if (r)=«f(r) + £ ™L £(A~%-(^.-^,). 

(12) 



Aji = Sji - M 



(13) 



M,- 



Pj( r )P»( r ) 
p(r) 



dr,i,j = 1, ...,m- 1. 



■ and v^f can be given as follows: 



C = / Pi(r)«£(r)dr. 



-j OCC p 

1=1 J 



X( r K*( r 'K( r >j( r ) 



r — r 



(14) 



(15) 



(16) 



Atomic calculations || show that in all situations, the 
KLI potential mimics the OEP results extremely well and 
even correctly preserves the property of integer discon- 
tinuity d]ll[]. Moreover, the KLI method is also easy 
to program. In this article we restrict ourself to various 
exchange-only methods, however, the extension of these 
methods to exchange-correlation case is straightforward 
using the best available correlation energy functional. 

For the constants Wi = v^f p — Vi (Eq. (9)) an al- 
ternative expression is proposed by Gritsenko et al. JToj ] 
in terms of orbital energies q. It follows from gauge 
invariance requirements, proper scaling and short range 
behaviour of the response part of Eq. (13) that Wj can 
only depend on an energy difference. In their work Eq. 
(13) is approximated by: 



,SSP 



(r) 



8^2' 



occ — 1 

3^" X! 

1=1 



^(r)| 2 



(17) 



We use the notation SSP (Slater's potential + step po- 
tential) for this exchange potential. The parameter 
is determined from the homogeneous electron gas (HEG) 
model and therefore Eq. ([l7]) is exact in the HEG limit 
|]l9| . This constant is chosen as universal parameter for 
all the calculations. The main advantage of this expres- 
sion is that one can avoid the matrix inversion of Eq. 
(12). Note that the summation runs over all the oc- 
cupied orbitals except the highest one like in Eq. (12). 
With the step potential-like second term v x SP provides a 
good approximation to the OEP exchange-potential [^9) . 

In actual calculations we used the following iterative 
scheme for getting v^ LI instead of Eq. (jl^) in each kth 
SCF cycle 



v KL I(k +, ){v) = v S ik){r) + P f k {v«"M _^). (18) 

i=l " 

The initial guess for Eq. ([t8]) is chosen as follows 



(19) 



We found this particular form of the initial guess is close 
to the KLI exchange potential. 
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The KLI approximation to OEP provides a self- 
interaction free exchange potential and proper — i 
asymptotics as r — > oo. The approximate potentials v~! 
and v% SP exhibit proper asymptotics as well . LDA or 
GGA exchange potentials do not show the correct asymp- 
totic behaviour JtJ. For instance, while the popular GGA 
exchange-correlation energy functional |2lf| provides en- 
ergies for atoms with nice agreement with the exact re- 
sults, the exchange-correlation potentials and the energy 
eigenvalues of the highest occupied state, are both suf- 
fering from error Nevertheless, various forms of 
GGAs are of great importance and further developments 
can certainly be expected in this field in the next future 
{§• 

3. Computational details 

Closed shell exchange-only KLI calculations have been 
performed for polyethylene and the results are compared 
with those from the Hartree-Fock and Slater approach. 
For all the calculations a new code is used which has been 
developed in our laboratory |33|. The various one and 
two-electronic integrals are taken from the Erlangen peri- 
odic Hartree-Fock program |}2j . During the calculations 
the number of fc-points is set to 25. The SCF iterations 
are continued until the density matrix elements changed 
by less then 10 -5 . 876 and 778 grid points are used for 
heavy atoms and for hydrogen in the numerical integra- 
tion using Becke's fuzzy cell method [Q. Thus the total 
number of grid points per unit cell is 2432. Further com- 
putational details of our implementation of polymeric cal- 
culations are given elsewhere [[l5|. For polyethylene the 
experimental "zig-zag" geometry is used Q. Since the 
various properties of polyethyene depend significantly on 
the geometry, we give the structure we have used for our 
calculations: the C-C and C-H bond lengths are set to 
1.54 and l.lOl while for the CCC, HCH and HCC an- 
gles 113, 108 and 113 degrees are used. Recent geometry 
optimization calculations with various DFT functionals 
provided values in very nice agreement with these ex- 
perimental bonding parameters (35| . Two types of bases 
are used: the dementi's minimal (7S/3P) and the larger 
double-C (9S/5P) basis sets for carbon and the 4S and 
6S for hydrogen atom, respectively The (11S/7P) 

basis set is also employed in certain cases. To make com- 
parison with previously published atomic and molecular 
results ||,[l2| and also to check the reliability of our code 
we have made calculations in the atomic and molecular 
limit. To do so, we set the translation vector to roughly 
20.0 a.u. and choose the interaction only between the 
nearest- neighbour CH2 unit cells Using these pa- 

rameters the computed properties must be very close 
to those obtained by atomic or molecular codes. The 
test provides nice agreement with atomic calculations ob- 



tained by Hartree-Fock, Slater's, KLI or SSP exchange- 
only method which confirms the reliability of our code. 
Further computational details of our implementation of 
KLI for polymers are given elsewhere [[l5| . 

4. Results 

The comparison is carried out for Be, LiH and for N2 
]l2| , p~9| | in the molecular limit and the energetical results 
are summarized in Table I. For the total energy we ob- 
tained -14.569, -14.562, -14.568 a.u. for Be using HF, 
Slater and KLI methods (with dementi's minimal ba- 
sis set). The deviation from HF can be compared with 
results obtained by others [PS: 7 and 1 mHartree for 
Slater and KLI, while Gritsenko et al. obtained 11 and 
0.0, respectively. For LiH the corresponding numbers are: 
-7.9543, -7.9542, -7.9542 a.u. The deviations: 0.0 , 0.0 
mH, while Grabo et al. [EJJ found 6 mH for Slater. The 
discrepancy between our and other results can be con- 
sidered rather small and is probably due to the different 
basis set applied here. For LiH the highest occupied or- 
bital energies are as follows: —.2914, —.3228 and —.3227 
a.u. which are compared with Grabo's values |0|, —.3017 
(HF), -.3150 (Slater) and -.3011 (KLI). Finally we give 
the test of our code for H2 as well. For two-electronic 
systems the results obtained by HF or KLI must be iden- 
tical ||i2|| . For Etot and e m (total energy and highest oc- 
cupied energy level) the values —1.12796,-1.127953 and 
-.5948, -.5913 a.u. are obtained by HF and KLI, re- 
spectively. The small difference must be attributed to 
numerical inaccuracy in the numerical integration (KLI). 
Results are also given for N2 in Table I. The KLI total 
energy for N2 differs from the HF one by 9 mH which 
is comparable with that obtained by a fully-numerical 
basis-set-free code ft2| . 

Our results for some physical properties of polyethy- 
lene are shown in Tables II-III. The results are listed in 
terms of the order of the neighbouring interaction. Actu- 
ally, it turns out that at least 5 neighbours are necessary 
in order to get convergent results for all the physical prop- 
erties we are interested in within the finite neighbouring 
approach jlq] . For comparison we have calculated all the 
properties at higher number of neighbouring cells |T5| ] 

Our calculation using the large (11S/7P) basis set 
provides a deeper Hartree-Fock total energy (—39.02326 
a.u.) than the one obtained by Suhai |44|] using the 6- 
31G* basis set (—39.02251 a.u.). According to the ex- 
pectations the comparison of total energies shows us that 
Hartree-Fock provides the deepest energies while the var- 
ious Kohn-Sham schemes result in somewhat higher en- 
ergies. In principle the total energy is always lower in 
HF than in (exchange-only) OEP since the use of a local 
potential imposes one more variational contraint on the 
orbitals B. In general, the following inequality holds: 
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E HF < Eoep < E KLI < E LDA @. KLI is always 
an upper bound to OEP and always lower than Elsdx 
(local exchange-only spin-density theory) . The difference 
is around 10 mH for KLI. For molecules T. Grabo and 
E. K. U. Gross [Q found differences around 7-8 mH for 
KLI. For a larger basis set KLI and HF total energies 
differ from each other more significantly (23 mH), which 
is much smaller, however, than the corresponding differ- 
ence for Slater (44 mH). We believe that the bulk part 
of the difference is not due to the approximate nature of 
the exchange potentials applied here, but to the different 
nature of Hartree-Fock and DFT approaches in accor- 
dance with earlier studies These results must be 
compared with the values obtained by calculations in the 
molecular limit |12|. Treating the CH<i unit cell in the 



molecular limit (Table I) , one can see that all the proper- 
ties (total energy, one electron energies) are close to the 
Hartree-Fock values. However, in polymer calculations 
we got for the one electron energies quite significant de- 
viation from the Hartree-Fock values (Tables II-III) : em- 
ployment of local exchange potentials generally leads to 
higher Fermi level and and at the same time to lower first 
virtual levels. Therefore, the single particle band gap of 
polyethylene is smaller for the KS methods than for HF. 

We also observe the peculiar behaviour of the highest 
occupied orbital energies e m obtained by KLI, which dif- 
fer from HF values significantly and are closer to [15| . 
This latter result is rather surprising since atomic and 
molecular calculations show that the e m values are very 
close to those obtained by HF or OEP Also, in 

HF we found no big difference between the highest occu- 
pied orbital energy in the molecular limit and in the infi- 
nite system. While this comes as a surprise we attribute 
the increase of the Fermi level due to periodic effects, 
which come into play only, when one employs localized 
exchange-potentials for infinite systems. One can follow 
how the periodic effect comes into play when the change 
of e m is examined in terms of the order of the neighbour- 
ing interaction in Table II-III. As may be read from Table 
II-III. , when neig = 1, e m obtained by Slater or KLI is 
the lowest one, however, when neig increased the Fermi 
level is getting closer to zero progressively. Hartree-Fock 
does not produce such a phenomenon since it has a non- 
local orbital dependent exchange-potential, which is al- 
ways deeper than the local counterpart, and thus is keep- 
ing the Fermi level at a deeper energy level. However, 
using the largest basis set in this study (11S/7P) we ob- 
tained somewhat higher e^ F (Table III) as well. As may 
be seen from Table III, by increasing the number of neigh- 
bours and the basis set e^ LI moves upward significantly 
and becomes higher than the corresponding LDA value 
of —0.242 a.u. except for SSP. As can be seen in 
Tables III, increasing the basis set reduces the calculated 
KS gap for all the methods, which is mainly due to the 
lowering of the first virtual level e m +i and to the upward 



movement of e m . While the HF single-particle gap is 
still too wide (16.6 eV), the calculated Slater's and KLI 
eigenvalue gaps are too narrow (3.4 and 5.1 eV). Others 
reported the value of 16.7 eV for the HF gap quite close 
to our value It is worth to note, however, that other 
calculations with 6-31G** basis set resulted in the smaller 
13.4 eV for the HF gap 0]. Recent reports proved that 
the band gap strongly depends not only on the basis set 
but also on the quality of the finite neigbour approach 
@- 

It is known (4^] that the exact exchange-correlation 
potential exhibits a discontinous jump as the number of 
electrons passes through an integer. This derivative dis- 
continuity plays an important role in the calculations of 
the exact band gap which can be written as 



A — A KS + A 

— '-^nonint ' ^xc 



(20) 



-N+ 



-N 



(N) 



where N is the number of electrons, A^?™"* is the Kohn- 
Sham band gap and A xc is the derivative discontinuity. 
For continuum approximations like LDA or GGA this dis- 
continuity vanishes. However, in OEP as well as in our 
(orbital-dependent) approximations to OEP, the discon- 
tinuity is finite ^|Jl5|,^2| . The exchange-only contribution 
to the band gap, A x , is given by P,[22| 



A x (i 



V) = (U v \V r , 



-HF 



{ Ui \v» F 



v° x EP [p]\u 

-OEP 



(21) 



v 



x 



u. 



{vi\vi) — 2{vi\iv), 



with v x F being the nonlocal HF exchange operator, how- 
ever, constructed from the N/2 occupied KS orbitals. 
In our previous study fl5|| we have calculated A x using 
Slater and SSP methods and the results can also be seen 
in Table II-III. The last two terms on the r.h.s. of Eq. 
(20) will vanish for systems with periodic boundary con- 
ditions and an infinite number of unit cells Since we 
only take a finite number of cells into account, we will 
also get contributions to A^. from these terms. We have 
studied the convergence of A x and, for a given basis set, 
at neig — 5 the discontinuity seems to be poorly con- 
verged with respect to the basis set and the results have 
to be read with due caution. In principle the discon- 
tinuity is guaranteed to be smaller than the true band 
gap by its definition However, calculations do not 
confirm this expectation for semiconductors Q. Never- 
theless, A x brings the KS gap somewhat closer to the HF 
one. With the minimal basis set A^ represents a rather 
small correction to the band gap and actually is negative 
(Table II). This is in contrast with the finding that A x 
amounts to typically twice the band gap ||. However, 
the situation changes completely when a larger basis set 
is used (Table III). The corrected band gap obtained by 
the larger basis set with the KLI method (Table III) is 
12.8 eV while the Slater and SSP approach yield 9.7 and 
12.7 eV, respectively 111. 
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The much more involved second-order self-energy cal- 
culations resulted in 10.3 eV for the quasi-particle (QP) 
band gap |p4 44|. Further corrections to A would cer- 
tainly result in a gap close to the QP and experimental 
values. According to the expectations KLI and SSP pro- 
vide values not far from the HF one while the naked 
Slater potential yields somewhat higher band gap then 
the experimental value of 8.8 eV |2^j3^|. In the basis 
set limit one expects the KLI gap to be close to the HF 
value. The discrepancy beetwen the KLI gap and the 
experimental one must be partly attributed to the corre- 
lation contribution to the derivative discontinuity (A c ). 
It is also useful to make comparison with quasi-particle 
calculations where one corrects the HF single-particle gap 
to get the quasi-particle band structure (2^^(],Q. Ac- 
cording to the results of Suhai obtained for polyethylene 
|Q the QP correction reduced the HF gap from 15.0 eV 
to 10.3 eV while the derivative discontinuity correction 
raises the OEP gap from 5.1 eV up to 12.8 eV. It would 
be rather interesting to add on further corrections to 
including correlation at second order level of theory [B2| . 
However, these kind of calculations are beyond the scope 
of the present article. 

Atomic calculations provided KLI eigenvalues which 
are very close to the HF and OEP values In 
molecules, however, Grabo et al. found slight deviations 
from the HF values jl2| . Approximate OEP calculations 
on polyethylene also resulted in eigenvalues different from 
HF 11511. In Table II-III we have already demonstrated 
that this is true in the KLI approach as well. In order to 
make a more detailed analysis of this problem we give the 
one-electron energies obtained by HF and KLI exchange- 
only methods in Table IV. We compare the HF and KLI 
eigenvalues obtained in the molecular limit and in the 
polymer. In general HF provides deeper eigenvalues by 
100 to 200 mH for the core orbitals and much less differ- 
ence for the valence orbitals. The reverse is true for the 
KLI virtual orbital energies. The KLI e m value of —0.093 
eV is compared with the values of —0.379 and —0.242 eV 
obtained by HF and LDA jlf|, respectively. To sum up, 
the KLI eigenvalues shift upward in the valence region 
and shift downward in the virtual space when compared 
with the HF values in the polymer. The KLI highest oc- 
cupied orbital energy is rather close to the HF one in the 
molecular limit. The differences are around 15 and 100 
mH, in the molecule and in the polymer. 

The exchange energies are also presented for various 
methods. KLI produces a higher exchange energy than 
HF by some 30 mH. The reason for this has to be found 
in the different form of the virial theorem which holds for 
HF and for DFT. In DFT the virial theorem is as follows 

E x [p] = J dr[3p(r)+rVp(r)K(r),. (22) 
On the basis of Eq. (G3) deeper exchange-potential will 



provide deeper exchange energy for the same density. 
Among the model potentials employed in this study the 
Slater potential exhibits the deepest potential curve in 
atoms 0,|ll[l and in Table III one can see that the ex- 
change energy E^ later is the deepest and differ from 
Hartree-Fock by some 50 mH. KLI provides exchange 
potential similar to that of LDA jl9],[llj in magnitude, 
although the shape of the curves are different (see e.g. 
in refs. [0jl|]). The Fermi level e m is deepest for the 
Slater's potential among the KS-based methods due to 
the deepest exchange potential. Even the first virtual 
level dropped below zero as obtained by the largest basis 
set by the Slater's exchange potential. The exchange en- 
ergy obtained by this potential is approximately by 1% 
lower than the Hartree-Fock exchange energy. This dif- 
ference is around the magnitude of the non-dynamical 
correlation energy JT(|. Therefore the Slater's potential 
is a suitable candidate for further calculations together 
with the commonly used correlation energy functionals 
which account for only the short range type dynami- 
cal correlation energy. The success of the popular ex- 
change energy functionals, such as the one of Becke |45] ] 
and the GGA's [|T) lie in the overestimation of Hartree- 
Fock exchange and therefore in the effective treatment 
of exchange-correlation mainly due to the cancellation 
of errors with opposite signs |Tq |. On the basis of Ta- 
ble III SP provides the best agreement either with the 
Hartree-Fock total energy or with the exchange energy. 
Therefore one might speculate that the SSP potential is 
somewhat closer to the exact OEP potential than to the 
KLI potential. The v^ SP is the functional not only of 
the orbitals but also of the eigenvalues of the KS prob- 
lem while v^ LI is only of the orbitals. Further steps in 
this direction can be desirable in the future in developing 
a general and accurate potential v^ p [{ui 7 e;}] ^ which 
is a functional of the eigensolutions of the KS problem. 

In Fig. 1 the band structure is plotted obtained by var- 
ious methods using the Clementi's double-C basis. One 
can see that virtual bands are appearing below the zero 
energy level obtained only by the Slater potential. Note 
that the first three curves in both of the figures are the 
occupied valence energy levels (Fig. 1 and 2, respectively, 
the lowest, core level is not plotted). We also note from 
Fig. 1 that the highest three virtual levels are similar to 
the virtuals obtained by the minimal basis set (see Fig. 
1 in ref. Jl5|). However, new patterns appear for the first 
two virtual levels as obtained by the larger basis set (plot- 
ted with dashed lines on Fig. 1). The lowest transition 
from the Fermi level to the first virtual is appearing at 
the edge of the Brillouin zone (fc = 1). However, around 
k = 0.5 one can see quasi-degeneracy of certain bands 
or even crossing of virtual levels close to the edge of the 
Brillouin zone. 

The maximum of the valence band is appearing in the 
range of [-2,-11] eV with the lowest and highest values 
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for KLI and Hartree-Fock, respectively, compared to the 
experimentally suggested values of ionization potential 
of 9.6 — 9.8 eV |38|. Others give lower experimental val- 
ues in the 7.6 — 8.8 eV range |^4|,^5|. KLI provides the 
much lower 2.3 eV value while Slater and SSP give 7.6 
and 4.4 eV jfj|, respectively. By density functional linear 
muffin-tin orbital method the value of 5.1 eV is obtained 
for helical polyethylene |25|. The calculated lowest va- 
lence band widths for Xa, HF, SSP, Slater's and KLI 
approach are 6.2, 9.0, 4.0, 3.6 and 3.9 eV respectively, 
which are to be compared with the experimental 7.2 eV 
Q . The total valence bandwidths are 14.0 eV {Xa), 
15.4 eV (SSP), 15.6 eV (KLI) and 19.8 eV (HF), com- 
pared to a experimental value of 16.2 eV |58| . Slater, KLI 
and SSP perform quite well for this band width. The 
calculated gap between the lowest valence band and the 
minimum of the higher valence bands is 3.3 eV (HF), 1.9 
eV (Xa), 3.0 (SSP), 3.2 (Slater) || and 2.8 eV (KLI) 
compared with the experimental 2.0 eV J38|. The HF 
and Xa values are in accordance with those obtained by 
others HSfl]. The bottom of the valence band with a 
symmetry is 18.0 eV obtained by KLI which is compared 
with the HF (31.0 eV) and Slater (23.6) eV values. The 
overall conclusion can be drawn that the various orbital 
dependent potentials provide band structures similar to 
one another and all of them differ from HF. The largest 
deviation from experiment is found for the lowest valence 
band width (LVBW) which is rather narrow in the OEP 
methods. HF give too wide LVBW while Xa-LDA yields 
a value in closest agreement with experiment. 



those obtained by Hartree-Fock. The addition of the 
step-potential as a response part of the exact exchange 
(Eq. 13) to the Slater's potential considerably improves 
the quality of various physical properties calculated for 
polyethylene. To test the quality and reliability of our 
code we performed calculations in the atomic and molec- 
ular limit. The results can be brought into harmony with 
those obtained by others. 

In general, we find that the band structure calculated 
by different v x are similar with little qualitative differ- 
ence, while all of them differ significantly from HF which 
have valence bands much deeper and much higher virtual 
levels. The single-particle band gap obtained by the KLI 
or Slater approaches is too narrow when compared with 
the experiment. The shifted values are observed for the 
Fermi level, which differ from Hartree-Fock values signifi- 
cantly. The one electron energy levels are carefully exam- 
ined and an upward shift of eigenvalues observed when 
compared with those of HF, and also the counter shift of 
virtual levels. The deviation from the HF eigenvalues is 
more significant then in molecules or in atoms. In order 
to make a valuable comparison with Hartree-Fock, LDA 
and with the experimental gap, correction must be taken 
into account to the Kohn-Sham eigenvalue gap. The ex- 
change contribution to the derivative discontinuity is cal- 
culated by the KLI exchange potential. In this way the 
KS single particle gap is corrected by a single shot calcu- 
lation on top of exchange-only OEP. The exchange contri- 
bution to the band gap is somewhat less than twice the 
KS eigenvalue gap. According to the expectations the 
corrected gap is much closer to the Hartree-Fock value 
than the KS eigenvalue gap. 



Conclusions 

Density functionals of the third generation treat both 
the kinetic and exchange energy exactly at the orbital de- 
pendent level of the theory. These orbitals come from a 
local potential and, due to the Hohenberg-Kohn theorem, 
are therefore functionals of the density. The optimized 
effective potential method, in the form given by Talman 
and Shadwick Q allows one to treat orbital-dependent 
functionals in the framework of Kohn-Sham density func- 
tional theory. The OEP integral equations are difficult to 
solve for extended systems and therefore several approx- 
imations are introduced. Among them the one obtained 
by Krieger et al. (KLI) and another one derived by Grit- 
senko et al. have become known in the last few years. 

In this work we have applied the exchange-only KLI 
approximation to the OEP approach for insulating poly- 
mers. Calculations have been carried out for polyethy- 
lene with a new polymer code using the efficient KLI 
approximation to OEP. The results obtained by the 
"naked" Slater's potential and KLI are compared with 
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TABLE I. Comparison of our HF and KLI results 
(Clementi's double zeta basis) using the code DIOGENES 
[B3J with values from the literature. 

Etot and e m are the total energy and the highest occupied 
orbital energy in a.u. 







jjj--pot/ier 


KU our 


KU other 


Be Etot 


-14.569 


-14.5730 a 


-14.5684 


-14.5723 b 


Cm 


-0.3089 


-0.3093 


-.3092 


-.3089 


LiH Etot 


-7.954 


-7.987 b 


-7.954 


-7.987 c 




-7.8629 d 


-7.8620 e 






Em 


-0.2914 


-0.3017 b 


-0.3227 


-0.3011 c 




-.2870 d 


-.2857 e 







N 2 


Etot 

fra - 


-108.892 ' 
-0.6352 f 


-108.994 d 
-0.6152 c 


-108.883 ' 
-0.6659 f 


-108.986 c 
-0.6818 c 


CH 2 


Etot 
£m 


-38.8650 ' 
-0.3750 f 




-38.8632 ' 
-0.3595 f 







d with STO-3G basis 

c with STO-3G basis using Gaussian94 program package 
The large 11S/7P basis set is used for N2 and for CH^. 



Fig. 1 The calculated valence band structure obtained by various 
exchange-only methods as a function of the dimcnsionlcss variable k 
with k — being the zone center and k — 1 the zone boundary. The 
Brillouin zone is that corresponding to a CH2 unit cell. Solid lines 
correspond to the occupied levels while the dashed lines to the virtual 
levels. The two lowest dashed curves of virtual levels arc of particular 
interest (see text). The Clementi's double-^ basis set is used and S 
neighbours are considered for the CH2 unit cell. 
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TABLE III. Calculated properties of polyethylene by vari- 
ous DFT methods using the dementi's double-zeta basis set 
(9S/5P). 

Results with basis (11S/7P) are also given in the bottom of 
the Table. The sum of the calculated KS eigenvalue gap 
^nonint and the exchange derivative discontinuity A x is also 
given together with the exchange-only LDA and experimental 
gap (eV). The notations are the same as in Table I. 



TABLE II. Calculated properties of polyethylene by vari- 
ous DFT methods using the minimal basis set 
HF, Slater, SSP and KLI denote the Hartree-Fock, ex- 
change-only method with Slater's potential, Slater's poten- 
tial with orbital dependent step potential and the KLI ex- 
change-only OEP method. All the properties are in a.u. 
except the HOMO-LUMO gap which is given in eV. The 
dementi's minimal basis set (7S/3P) is used. Neig gives the 
number of neighbouring cells taken into account. Etot, E x are 
the calculated total and exchange energies per CH2 unit cell. 
A x is the exchange component of the derivative discontinuity. 
The calculations are carried out for the CH2 unit cell. 



neig 


HF 


Slater a 


SSP a 


KLI 


1 Etot (a.u.) 


-39.03535 


-39.01888 


-39.02768 


-39.02042 


E x (a.u.) 


-5.82128 


-5.84774 


-5.79080 


-5.7713 


e m (a.u.) 


-.6672 


-.5634 


-.4924 


-.3524 


em+i (a.u.) 


.3978 


.0074 


-.1084 


.2228 


gap (eV) 


29.0 


15.1 


16.4 


15.7 


A x (eV) 




-5.7 


-4.4 


-3.2 


5 


-38.88012 


-38.86891 


-38.87278 


-38.86981 




-5.75048 


-5.76397 


-5.7379 


-5.7297 




-.4747 


-.3524 


-.2312 


-.0871 




.3444 


-.0433 


.1310 


.2786 




22.3 


8.4 


9.9 


10.0 






-3.5 


-1.9 


-1.2 


8 


-38.88008 


-38.86890 


-38.87275 


-38.86954 




-5.75049 


-5.76400 


-5.7378 


-5.7284 




-.4760 


-.3518 


-.2318 


-.0972 




.3436 


-.0425 


.1303 


.275 




22.3 


8.4 


9.9 


10.1 






-3.5 


-1.9 


-1.1 


KS gap+Az 




4.9 


8.0 


9.0 



neig 


HF 


Slater a 


SSP a 


KLI 


5 Etot 


-39.01068 


-38.96683 


-38.99704 


-38.99311 


Ej; 


-5.8783 


-5.9368 


-5.8841 


-5.8405 


£771 


-.4049 


-.2951 


-.1628 


-.1023 


£m+l 


.1489 


-.1590 


-.0007 


.0284 


KS gap 


15.1 


3.7 


4.4 


3.6 






3.4 


3.8 


1.9 


8 (9S/5P) 


-39.01056 


-38.96684 


-38.99760 


-38.99407 




-5.8780 


-5.9367 


-5.9036 


-5.8621 




-.3984 


-.2774 


-.1130 


-.0997 




.1558 


-.1655 


-.0755 


.0322 




15.1 


3.7 


4.4 


3.6 






3.2 


3.7 


2.4 


8 (11S/7P) 


-39.02326 


-38.99926 


-39.01361 


-39.00701 




-5.8977 


-5.9486 


-5.9036 


-5.8553 




-.3786 


-.2401 


-.1130 


-.0933 




.2320 


-.1151 


.0755 


.0958 




16.6 


3.4 


5.1 


5.1 






6.3 


7.6 


7.7 


KS gap+Az 




9.7 


12.7 


12.8 


LDA 








7.6 a 


EXP 








8.8 b 



a '-l-l 

b |2J3ij] 



TABLE IV. One-electron energies (a.u.) obtained by 
HF and KLI exchange-only approaches using the dementi's 
(11S/7P) basis set and 8 neighbours for polyethylene (k — 1) 







KLI moI 




KLI poly 


£1 


-11.289 


-10.201 


-11.219 


-9.820 


£2 


-.904 


-.774 


-.800 


-.511 


£3 


-.587 


-.537 


-.681 


-.401 


£4 


-.375 


-.359 


-.379 


-.093 


£5 


.056 


-.318 


.232 


.096 


£6 


.240 


-.058 


.232 


.109 


£7 


.314 


-.005 


.392 


.301 


£8 


.518 


.192 


.719 


.527 
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The optimized effective potential (OEP) method allows for 
calculation of the local, effective single particle potential of 
density functional theory for explicitly orbital-dependent ap- 
proximations to the exchange-correlation energy functional. 
In the present work the OEP method is used together with 
the approximation due to Krieger, Li and Iafrate (KLI). We 
present the first application of this method to polymers. KLI 
calculations have been performed for the insulating polyethy- 
lene and the results have been compared to those from other 
orbital-dependent potentials. Various properties of the band 
structure are also calculated. The single-particle band gap 
strongly depends on the basis set with larger basis sets yield- 
ing narrow gaps. For certain physical quantities such as the 
total energy and the exchange energy, the various orbital- 
dependent exchange-only and Hartree-Fock results differ only 
slightly. For the highest occupied orbital energy the differ- 
ence is more pronounced. In order to get the right band gap 
in OEP the exchange contribution to the derivative discon- 
tinuity is calculated and added to the Kohn-Sham gap. The 
corrected gap obtained by the KLI approach is 12.8 eV com- 
pared with the Hartree-Fock and experimental values of 16.6 
and 8.8 eV, respectively. We observe, however, the strong 
dependence of the derivative discontinuity on the basis set. 



1. Introduction 

The optimized effective potential (OEP) method was 
invented by Talman and Shadwick [1], following the orig- 
inal idea of Sharp and Horton [2]. In the now classic 
OEP approach a given total energy functional depend- 
ing explicitly on single-particle orbitals is minimized un- 
der the constraint that these orbitals are solutions of a 
single-particle Schrodinger equation with a local, effective 
potential. The exchange-only OEP can be interpreted 
as the implementation of the exact exchange-only den- 
sity funtional theory (DFT) [3,4]. Very recently OEP is 
extended to include electron correlation effects via per- 
turbation theory [5,6]. In this theory both the correla- 
tion energy and the potential are expanded using stan- 
dard perturbation theory. In this way a self-consistent 
parameter-free KS ab initio scheme can be given. The 
exact energy functional and the exchange-correlation po- 
tential provided by the Hohenberg-Kohn theorem can be 
calculated order-by-order by means of the OEP pertur- 
bation theory [6]. 

To the best of our knowledge, full OEP calculations 



have so far only been performed for atoms [4,8] and for 
semiconductors [9] at the exchange-only level of the the- 
ory. The incorporation of electron correlation effects into 
the OEP formalism has been the subject of few studies 
[4,10]. Most of these studies, however, have used the 
complicated self-energy formalism. 

To reduce the numerical complexity of the OEP 
method, Krieger, Li and Iafrate (KLI) introduced a re- 
markably accurate analytical approximation to OEP and 
transformed the original OEP equations into an easily 
manageable form [7] . Applications to atoms and to solids 
show that the KLI method yields results that are nearly 
identical to those of OEP [8,11]. The KLI method has 
been successfully applied for molecules [12] and for cer- 
tain semiconductors [13,14] using pseudopotentials. 

In our previous study we have carried out calculations 
by means of various orbital-dependent exchange poten- 
tials which are accurate approximations of OEP and re- 
quire less computational effort [15]. We found that the 
various orbital dependent potentials provide a narrow 
DFT gap for polyethylene (PE). KS theory will typically 
give too narrow valence bands and too compressed con- 
duction bands [11]. The KS theory will only give an exact 
result for the highest occupied energy eigenvalue, there- 
fore, calculated states will in general not have the desired 
accuracy as one might hope [3,7]. In addition to this, it 
has been shown that the difference between the highest 
occupied and lowest unoccupied DFT eigenvalues in the 
iV-electron system (which is what is usually called the 
DFT band gap) is not the true quasiparticle band gap, 
but differs from it by the discontinuity in the exchange- 
correlation potential when an electron is added to the 
system [4,23]. Therefore, the so-called derivative discon- 
tinuity corrections must be taken into account in order to 
get the right gap [4,15,22]. Consequently, it is desirable 
to calculate the band structure using the KS method, and 
then simply add on a single-shot correction to all states 
[11,15]. In this way much better agreement is found with 
the Hartree-Fock gap when the exchange corrections are 
considered to the KS eigenvalue band gap [15]. None of 
the state-of-the-art exchange-correlation functionals de- 
pending only explicitly on the density show the above- 
mentioned discontinuity in the exchange-correlation po- 
tential and therefore do not yield the correct values for 
the band gaps of the insulating solids [11]. 

In the last decades quite a number of publications 
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appeared on polymer electronic structure using mainly 
Hartree-Fock ab initio methods [24]. Results are reported 
only in recent years, however, using DFT methods like 
local density approximation (LDA) [26]. The general- 
ized gradient approximation (GGA) [21] resulted in only 
modest improvement over LDA gap in semiconductors 
[27]. 

In this paper, we briefly review the theory of OEP and 
KLI methods. This work can be taken as the systematic 
continuation of our previous work on polymer electronic 
structure obtained by various orbital dependent local ex- 
change potentials [15]. Particular attention is paid to 
the correction of the OEP eigenvalue band gap. We will 
present the first results obtained by the KLI orbital de- 
pendent exchange potential for polymers. We give the 
calculated band gaps, band widths, band structures and 
total energy of polyethylene within the Kohn-Sham den- 
sity functional scheme using various orbital dependent 
exchange potentials. 



2. Basic Formalism 

The following spin-restricted Kohn-Sham type mean- 
field approach is employed for polymers (the extension of 
the theory for spinpolarized problems is straightforward) : 



[-^V 2 + v s (r)]ui(r) = ejUi(r) 



(1) 



where {uj(r)} are the single particle orbitals and v s (r) is 
the Kohn-Sham effective single particle potential [3]. 



v s (r) = v ext (v) + v H (r) + v xc (r) 



(2) 



The Kohn-Sham one-body potential contains the exter- 
nal, Hartree and exchange-correlation potential. v ex t(r) 
represents the Coulomb potential of the nuclei. vn(r) is 
the classical Coulomb repulsion of the electron charge 
density.The exchange-correlation potential v xc (r) is a 
nonclassical term. It is formally defined as functional 
derivative of the exchange-correlation energy E xc [p\: 



= (r,[p]) 



SE xc [p] 
Sp(r) 



(3) 



We seek an expression which contains the exact exchange 
potential and correlation can be added easily on top of 
the exact exchange-only theory using approximate cor- 
relation energy functionals [3]. It must be emphasized, 
however that the correlation energy functionals either in 
their local or gradient corrected form may not work in 
combination with exact exchange. This is because basi- 
cally all common correlation functionals do not have a 
long-range component in the corresponding correlation 
hole (the combined exchange-correlation hole is typically 
short-ranged) . The Colle-Salvetti gradient corrected cor- 
relation functional [28], which provided excellent results 



for atoms, performed rather badly for molecules for the 
abovementioned reason [8,12]. Approximate correlation 
functionals which are derived from the homogeneous 
or inhomogeneous electron gas model or obtained non- 
empirically from sum-rule conditions [16] have an incor- 
rect long-range tail (only dynamical correlation is ac- 
counted for). The correlation energy functionals provide 
improper local behaviour in the bond midpoint regions 
in molecules [16,19,20] which is mainly due to the lack 
of non-dynamical correlation effects [16]. Therefore, in 
those methods which work with exact exchange (HF and 
OEP) the correlation effect can only be accounted for 
in a more complicated way than in the so-called second 
generation of DFT [15]. 

In the optimized effective potential method (OEP) in 
contrast to the ordinary DFT, the exchange-correlation 
energy is approximated as an explicit functional of or- 
bitals and only implicit functional of the electronic den- 
sity [12]. The starting point of the OEP method is the 
total energy functional 

C P [P] = E/*<W(-^ 2 )^) 
+ j dvp(v)v ext {v) + \ J drdr'^^l+Eg^iiu,}]. (4) 

We make use of the chain rule for functional derivatives 
to obtain from Eq. (3) 



, 6E^[{ Ui }] 6 Ui (r>) 
Sui(r) Sp(r) 



+ c.c. (5) 



i=l J 



where formally are the orbitals which are, how- 

ever, implicit functionals of the density. Applying the 
functional chain rule again and after some algebra one 
obtains the following integral equation: 

occ ~ 

5>i(r) J dr'[v xc (r') - Vi (r')]G Bi (r, r'K(r') + c.c. = 



where 



and 



G si (r,r') = J2 



1 SE° c EP [ Ui ] 
uj(r) Uj(r) 



e« - e& 



(6) 



(7) 



(8) 



The integral in Eq. (6) is the fundamental expression 
for v xc in OEP and can be shown to be equivalent 
to Eq. (3) [4]. There is no known analytic solution 
for [{««}]• Therefore, only numerical solutions are 
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available for spherical atoms [1,4,7,29,30] and for solids 
[9,13,14] using atomiclike wave functions and assuming 
integer atomic occupation numbers within each muffin- 
tin sphere (linear augmented plane- wave method) [11]. 
These numerical solutions of Eq. (6) are confined to the 
exchange-only OEP and correlation has been taken into 
account only via approximate local functionals [9]. In 
the work of Krieger and co-workers the OEP integral 
equation is analyzed and a simple approximation is made 
which reduces the complexity of the original OEP equa- 
tion significantly and at the same time keeps many of 
the essential properties of OEP unchanged [7]. Krieger, 
Li and Iafrate gave an exact expression transforming the 
OEP integral equation (6) into a manageable form. They 
obtained the following, still exact expression for v xc (r): 



occ—l / \ 

C EP (r)=v S xc (r)+ £ E4EltfBP- Vi) + 



(9) 



1 occ 



V[^(r)V^(r)] 
p(r) 



J drpi(r)vi(r), 



where m = occ is the highest occupied one electron en- 
ergy level, pi = \ui\ 2 is the partial density. v s (r) is 
the Slater's potential (its correlation part is not known 
exactly) 
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2p(r) 



fir 



,l7(r,r') 



(10) 



where 7(1-, r ) is the first-order density matrix. The sum- 
mation runs over the orbital index for all the occupied 
orbitals up to the highest occupied mth orbital (Fermi 
level). The function pi is defined by 
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Ui(r) 



dr'l 



,OEP 



(r')-^(r')]G si (r,r'K(r'). 

(11) 



In practical applications the last term in Eq. (9) turned 
out to be quite small in atomic systems and has small 
effect only on the atomic shell boundaries [7]. This ad- 
ditional term's average over the mth orbitals is zero [7]. 
Neglecting this term one obtains the exchange-only KLI- 
approximation which, after some algebra, can be written 
in the following form [7] 



/ 



dr,i,j = 1, ...,m- 1. 



V s . a and v^F can be given as follows: 



-. occ ~ 

£/ 

i=l J 



dvd\ 



,U*{t)u*(t )Ui(v')Uj{T) 



(14) 



(15) 



(16) 



Atomic calculations [8] show that in all situations, the 
KLI potential mimics the OEP results extremely well and 
even correctly preserves the property of integer discon- 
tinuity [4,11]. Moreover, the KLI method is also easy 
to program. In this article we restrict ourself to various 
exchange-only methods, however, the extension of these 
methods to exchange-correlation case is straightforward 
using the best available correlation energy functional. 

For the constants w, = v®® 1 * — Vi (Eq. (9)) an al- 
ternative expression is proposed by Gritsenko et al. [19] 
in terms of orbital energies e^. It follows from gauge 
invariance requirements, proper scaling and short range 
behaviour of the response part of Eq. (13) that Wj can 
only depend on an energy difference. In their work Eq. 
(13) is approximated by: 



,SSP 



K(r)| 
p(v) 



(17) 



We use the notation SSP (Slater's potential + step po- 
tential) for this exchange potential. The parameter 
is determined from the homogeneous electron gas (HEG) 
model and therefore Eq. (17) is exact in the HEG limit 
[19]. This constant is chosen as universal parameter for 
all the calculations. The main advantage of this expres- 
sion is that one can avoid the matrix inversion of Eq. 
(12). Note that the summation runs over all the oc- 
cupied orbitals except the highest one like in Eq. (12). 
With the step potential-like second term v s SP provides a 
good approximation to the OEP exchange-potential [19]. 

In actual calculations we used the following iterative 
scheme for getting v x LI instead of Eq. (12) in each kth 
SCF cycle 



m — 1 



i—1 



.KLI 



(r) 



P(r) 



|2 occ—l 

L £ (a-%c 

J'=l 



Sji 



M 4 



(12) 



(13) 



The initial guess for Eq. (18) is chosen as follows 
«f" (0) (r) = |(t,^(r) + |t,f(r)). 



(19) 



We found this particular form of the initial guess is close 
to the KLI exchange potential. 



3 



The KLI approximation to OEP provides a self- 
interaction free exchange potential and proper — ^ 
asymptotics as r — > oo. The approximate potentials 
and v^ SP exhibit proper asymptotics as well [19]. LDA or 
GGA exchange potentials do not show the correct asymp- 
totic behaviour [7]. For instance, while the popular GGA 
exchange-correlation energy functional [21] provides en- 
ergies for atoms with nice agreement with the exact re- 
sults, the exchange-correlation potentials and the energy 
eigenvalues of the highest occupied state, are both suf- 
fering from error [6,7]. Nevertheless, various forms of 
GGAs are of great importance and further developments 
can certainly be expected in this field in the next future 
[39]. 

3. Computational details 

Closed shell exchange-only KLI calculations have been 
performed for polyethylene and the results are compared 
with those from the Hartree-Fock and Slater approach. 
For all the calculations a new code is used which has been 
developed in our laboratory [33]. The various one and 
two-electronic integrals are taken from the Erlangen peri- 
odic Hartree-Fock program [32]. During the calculations 
the number of A;-points is set to 25. The SCF iterations 
are continued until the density matrix elements changed 
by less then 10~ 5 . 876 and 778 grid points are used for 
heavy atoms and for hydrogen in the numerical integra- 
tion using Becke's fuzzy cell method [45]. Thus the total 
number of grid points per unit cell is 2432. Further com- 
putational details of our implementation of polymeric cal- 
culations are given elsewhere [15]. For polyethylene the 
experimental "zig-zag" geometry is used [26]. Since the 
various properties of polyethyene depend significantly on 
the geometry, we give the structure we have used for our 
calculations: the C-C and C-H bond lengths are set to 
1.54 and l.lOA while for the CCC, HCH and HCC an- 
gles 113, 108 and 113 degrees are used. Recent geometry 
optimization calculations with various DFT functionals 
provided values in very nice agreement with these ex- 
perimental bonding parameters [35]. Two types of bases 
are used: the Clementi's minimal (7S/3P) and the larger 
double-^ (9S/5P) basis sets for carbon and the 4S and 
6S for hydrogen atom, respectively [36]. The (11S/7P) 
basis set is also employed in certain cases. To make com- 
parison with previously published atomic and molecular 
results [8,12] and also to check the reliability of our code 
we have made calculations in the atomic and molecular 
limit. To do so, we set the translation vector to roughly 
20.0 a.u. and choose the interaction only between the 
nearest-neighbour CH2 unit cells [24]. Using these pa- 
rameters the computed properties must be very close 
to those obtained by atomic or molecular codes. The 
test provides nice agreement with atomic calculations ob- 



tained by Hartree-Fock, Slater's, KLI or SSP exchange- 
only method which confirms the reliability of our code. 
Further computational details of our implementation of 
KLI for polymers are given elsewhere [15]. 

4. Results 

The comparison is carried out for Be, LiH and for 
[12,19] in the molecular limit and the energetical results 
are summarized in Table I. For the total energy we ob- 
tained -14.569, -14.562, -14.568 a.u. for Be using HF, 
Slater and KLI methods (with Clementi's minimal ba- 
sis set). The deviation from HF can be compared with 
results obtained by others [19]: 7 and 1 mHartree for 
Slater and KLI, while Gritsenko et al. obtained 11 and 
0.0, respectively. For LiH the corresponding numbers are: 
-7.9543, -7.9542, -7.9542 a.u. The deviations: 0.0 , 0.0 
mH, while Grabo et al. [12] found 6 mH for Slater. The 
discrepancy between our and other results can be con- 
sidered rather small and is probably due to the different 
basis set applied here. For LiH the highest occupied or- 
bital energies are as follows: —.2914, —.3228 and —.3227 
a.u. which are compared with Grabo's values [12], —.3017 
(HF), -.3150 (Slater) and -.3011 (KLI). Finally we give 
the test of our code for H2 as well. For two-electronic 
systems the results obtained by HF or KLI must be iden- 
tical [12]. For E tot and e m (total energy and highest oc- 
cupied energy level) the values —1.12796,-1.127953 and 
-.5948, -.5913 a.u. are obtained by HF and KLI, re- 
spectively. The small difference must be attributed to 
numerical inaccuracy in the numerical integration (KLI) . 
Results are also given for N 2 in Table I. The KLI total 
energy for N2 differs from the HF one by 9 mH which 
is comparable with that obtained by a fully-numerical 
basis-set-free code [12]. 

Our results for some physical properties of polyethy- 
lene are shown in Tables II-III. The results are listed in 
terms of the order of the neighbouring interaction. Actu- 
ally, it turns out that at least 5 neighbours are necessary 
in order to get convergent results for all the physical prop- 
erties we are interested in within the finite neighbouring 
approach [15]. For comparison we have calculated all the 
properties at higher number of neighbouring cells [15] 

Our calculation using the large (11S/7P) basis set 
provides a deeper Hartree-Fock total energy (—39.02326 
a.u.) than the one obtained by Suhai [44] using the 6- 
31G* basis set (-39.02251 a.u.). According to the ex- 
pectations the comparison of total energies shows us that 
Hartree-Fock provides the deepest energies while the var- 
ious Kohn-Sham schemes result in somewhat higher en- 
ergies. In principle the total energy is always lower in 
HF than in (exchange-only) OEP since the use of a local 
potential imposes one more variational contraint on the 
orbitals [4]. In general, the following inequality holds: 
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E HF < Eoep < E KLI < E LDA [18]. KLI is always 
an upper bound to OEP and always lower than Elsdx 
(local exchange-only spin-density theory) . The difference 
is around 10 mH for KLI. For molecules T. Grabo and 
E. K. U. Gross [12] found differences around 7-8 mH for 
KLI. For a larger basis set KLI and HF total energies 
differ from each other more significantly (23 mH), which 
is much smaller, however, than the corresponding differ- 
ence for Slater (44 mH). We believe that the bulk part 
of the difference is not due to the approximate nature of 
the exchange potentials applied here, but to the different 
nature of Hartree-Fock and DFT approaches in accor- 
dance with earlier studies [7,12]. These results must be 
compared with the values obtained by calculations in the 
molecular limit [12]. Treating the CHi unit cell in the 
molecular limit (Table I) , one can see that all the proper- 
ties (total energy, one electron energies) are close to the 
Hartree-Fock values. However, in polymer calculations 
we got for the one electron energies quite significant de- 
viation from the Hartree-Fock values (Tables II-III) : em- 
ployment of local exchange potentials generally leads to 
higher Fermi level and and at the same time to lower first 
virtual levels. Therefore, the single particle band gap of 
polyethylene is smaller for the KS methods than for HF. 

We also observe the peculiar behaviour of the highest 
occupied orbital energies e TO obtained by KLI, which dif- 
fer from HF values significantly and are closer to e^ a [15]. 
This latter result is rather surprising since atomic and 
molecular calculations show that the e TO values are very 
close to those obtained by HF or OEP [7,8,12]. Also, in 
HF we found no big difference between the highest occu- 
pied orbital energy in the molecular limit and in the infi- 
nite system. While this comes as a surprise we attribute 
the increase of the Fermi level due to periodic effects, 
which come into play only, when one employs localized 
exchange-potentials for infinite systems. One can follow 
how the periodic effect comes into play when the change 
of e m is examined in terms of the order of the neighbour- 
ing interaction in Table II-III. As may be read from Table 
II-III., when neig = 1, e m obtained by Slater or KLI is 
the lowest one, however, when neig increased the Fermi 
level is getting closer to zero progressively. Hartree-Fock 
does not produce such a phenomenon since it has a non- 
local orbital dependent exchange-potential, which is al- 
ways deeper than the local counterpart, and thus is keep- 
ing the Fermi level at a deeper energy level. However, 
using the largest basis set in this study (11S/7P) we ob- 
tained somewhat higher e^ F (Table III) as well. As may 
be seen from Table III, by increasing the number of neigh- 
bours and the basis set e^ LI moves upward significantly 
and becomes higher than the corresponding LDA value 
of —0.242 a.u. [15], except for SSP. As can be seen in 
Tables III, increasing the basis set reduces the calculated 
KS gap for all the methods, which is mainly due to the 
lowering of the first virtual level e m+ i and to the upward 



movement of e TO . While the HF single-particle gap is 
still too wide (16.6 eV), the calculated Slater's and KLI 
eigenvalue gaps are too narrow (3.4 and 5.1 eV). Others 
reported the value of 16.7 eV for the HF gap quite close 
to our value [40]. It is worth to note, however, that other 
calculations with 6-31G** basis set resulted in the smaller 
13.4 eV for the HF gap [44]. Recent reports proved that 
the band gap strongly depends not only on the basis set 
but also on the quality of the finite neigbour approach 
[40]. 

It is known [42] that the exact exchange-correlation 
potential exhibits a discontinous jump as the number of 
electrons passes through an integer. This derivative dis- 
continuity plays an important role in the calculations of 
the exact band gap which can be written as 

A = A« s nint + A, c (20) 
= e« s +1 (N)-e« s (N)+A xc , 

where TV is the number of electrons, A^?g mi is the Kohn- 
Sham band gap and A xc is the derivative discontinuity. 
For continuum approximations like LDA or GGA this dis- 
continuity vanishes. However, in OEP as well as in our 
(orbital-dependent) approximations to OEP, the discon- 
tinuity is finite [4,15,22]. The exchange-only contribution 
to the band gap, A x , is given by [9,22] 

A x (i ->«,) = (u v \v? F - v° EP \p]\u v ) (21) 

— {vi\vi) — 2{vi\iv), 

with v x F being the nonlocal HF exchange operator, how- 
ever, constructed from the N/2 occupied KS orbitals. 
In our previous study [15] we have calculated using 
Slater and SSP methods and the results can also be seen 
in Table II-III. The last two terms on the r.h.s. of Eq. 
(20) will vanish for systems with periodic boundary con- 
ditions and an infinite number of unit cells [9] . Since we 
only take a finite number of cells into account, we will 
also get contributions to A^ from these terms. We have 
studied the convergence of A^ and, for a given basis set, 
at neig = 5 the discontinuity seems to be poorly con- 
verged with respect to the basis set and the results have 
to be read with due caution. In principle the discon- 
tinuity is guaranteed to be smaller than the true band 
gap by its definition [4]. However, calculations do not 
confirm this expectation for semiconductors [9]. Never- 
theless, A x brings the KS gap somewhat closer to the HF 
one. With the minimal basis set A^ represents a rather 
small correction to the band gap and actually is negative 
(Table II). This is in contrast with the finding that A^ 
amounts to typically twice the band gap [9]. However, 
the situation changes completely when a larger basis set 
is used (Table III). The corrected band gap obtained by 
the larger basis set with the KLI method (Table III) is 
12.8 eV while the Slater and SSP approach yield 9.7 and 
12.7 eV, respectively [15]. 
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The much more involved second-order self-energy cal- 
culations resulted in 10.3 eV for the quasi-particle (QP) 
band gap [24,44]. Further corrections to A would cer- 
tainly result in a gap close to the QP and experimental 
values. According to the expectations KLI and SSP pro- 
vide values not far from the HF one while the naked 
Slater potential yields somewhat higher band gap then 
the experimental value of 8.8 eV [24,38]. In the basis 
set limit one expects the KLI gap to be close to the HF 
value. The discrepancy beetwen the KLI gap and the 
experimental one must be partly attributed to the corre- 
lation contribution to the derivative discontinuity (A c ). 
It is also useful to make comparison with quasi-particle 
calculations where one corrects the HF single-particle gap 
to get the quasi-particle band structure [24,40,44]. Ac- 
cording to the results of Suhai obtained for polyethylene 
[44] the QP correction reduced the HF gap from 15.0 eV 
to 10.3 eV while the derivative discontinuity correction 
raises the OEP gap from 5.1 eV up to 12.8 eV. It would 
be rather interesting to add on further corrections to A x 
including correlation at second order level of theory [22]. 
However, these kind of calculations are beyond the scope 
of the present article. 

Atomic calculations provided KLI eigenvalues which 
are very close to the HF and OEP values [8]. In 
molecules, however, Grabo et al. found slight deviations 
from the HF values [12]. Approximate OEP calculations 
on polyethylene also resulted in eigenvalues different from 
HF [15]. In Table II-III we have already demonstrated 
that this is true in the KLI approach as well. In order to 
make a more detailed analysis of this problem we give the 
one-electron energies obtained by HF and KLI exchange- 
only methods in Table IV. We compare the HF and KLI 
eigenvalues obtained in the molecular limit and in the 
polymer. In general HF provides deeper eigenvalues by 
100 to 200 mH for the core orbitals and much less differ- 
ence for the valence orbitals. The reverse is true for the 
KLI virtual orbital energies. The KLI e m value of —0.093 
eV is compared with the values of —0.379 and —0.242 eV 
obtained by HF and LDA [15], respectively. To sum up, 
the KLI eigenvalues shift upward in the valence region 
and shift downward in the virtual space when compared 
with the HF values in the polymer. The KLI highest oc- 
cupied orbital energy is rather close to the HF one in the 
molecular limit. The differences are around 15 and 100 
mH, in the molecule and in the polymer. 

The exchange energies are also presented for various 
methods. KLI produces a higher exchange energy than 
HF by some 30 mH. The reason for this has to be found 
in the different form of the virial theorem which holds for 
HF and for DFT. In DFT the virial theorem is as follows 
[37], 

E x [p] = J dv[3p(r) + rVp(r)K(r), . (22) 
On the basis of Eq. (22) deeper exchange-potential will 



provide deeper exchange energy for the same density. 
Among the model potentials employed in this study the 
Slater potential exhibits the deepest potential curve in 
atoms [7,11] and in Table III one can see that the ex- 
change energy E x later is the deepest and differ from 
Hartree-Fock by some 50 mH. KLI provides exchange 
potential similar to that of LDA [19,11] in magnitude, 
although the shape of the curves are different (see e.g. 
in refs. [7,11]). The Fermi level e TO is deepest for the 
Slater's potential among the KS-based methods due to 
the deepest exchange potential. Even the first virtual 
level dropped below zero as obtained by the largest basis 
set by the Slater's exchange potential. The exchange en- 
ergy obtained by this potential is approximately by 1% 
lower than the Hartree-Fock exchange energy. This dif- 
ference is around the magnitude of the non-dynamical 
correlation energy [16]. Therefore the Slater's potential 
is a suitable candidate for further calculations together 
with the commonly used correlation energy functionals 
which account for only the short range type dynami- 
cal correlation energy. The success of the popular ex- 
change energy functionals, such as the one of Becke [45] 
and the GGA's [21] lie in the overestimation of Hartree- 
Fock exchange and therefore in the effective treatment 
of exchange-correlation mainly due to the cancellation 
of errors with opposite signs [16]. On the basis of Ta- 
ble III i>f SP provides the best agreement either with the 
Hartree-Fock total energy or with the exchange energy. 
Therefore one might speculate that the SSP potential is 
somewhat closer to the exact OEP potential than to the 
KLI potential. The v x sp is the functional not only of 
the orbitals but also of the eigenvalues of the KS prob- 
lem while v x LI is only of the orbitals. Further steps in 
this direction can be desirable in the future in developing 
a general and accurate potential v°f p [{ui, e,}] [6] which 
is a functional of the eigensolutions of the KS problem. 

In Fig. 1 the band structure is plotted obtained by var- 
ious methods using the Clementi's double-^ basis. One 
can see that virtual bands are appearing below the zero 
energy level obtained only by the Slater potential. Note 
that the first three curves in both of the figures are the 
occupied valence energy levels (Fig. 1 and 2, respectively, 
the lowest, core level is not plotted). We also note from 
Fig. 1 that the highest three virtual levels are similar to 
the virtuals obtained by the minimal basis set (see Fig. 
1 in ref. [15]). However, new patterns appear for the first 
two virtual levels as obtained by the larger basis set (plot- 
ted with dashed lines on Fig. 1). The lowest transition 
from the Fermi level to the first virtual is appearing at 
the edge of the Brillouin zone (k = 1). However, around 
k = 0.5 one can see quasi-degeneracy of certain bands 
or even crossing of virtual levels close to the edge of the 
Brillouin zone. 

The maximum of the valence band is appearing in the 
range of [-2,-11] eV with the lowest and highest values 



6 



for KLI and Hartree-Fock, respectively, compared to the 
experimentally suggested values of ionization potential 
of 9.6 — 9.8 eV [38]. Others give lower experimental val- 
ues in the 7.6 — 8.8 eV range [24,25]. KLI provides the 
much lower 2.3 eV value while Slater and SSP give 7.6 
and 4.4 eV [15], respectively. By density functional linear 
muffin-tin orbital method the value of 5.1 eV is obtained 
for helical polyethylene [25]. The calculated lowest va- 
lence band widths for Xa, HF, SSP, Slater's and KLI 
approach are 6.2, 9.0, 4.0, 3.6 and 3.9 eV respectively, 
which are to be compared with the experimental 7.2 eV 
[38]. The total valence bandwidths are 14.0 eV (Xa), 
15.4 eV (SSP), 15.6 eV (KLI) and 19.8 eV (HF), com- 
pared to a experimental value of 16.2 eV [38]. Slater, KLI 
and SSP perform quite well for this band width. The 
calculated gap between the lowest valence band and the 
minimum of the higher valence bands is 3.3 eV (HF), 1.9 
eV (Xa), 3.0 (SSP), 3.2 (Slater) [15] and 2.8 eV (KLI) 
compared with the experimental 2.0 eV [38]. The HF 
and Xa values are in accordance with those obtained by 
others [24,26]. The bottom of the valence band with a 
symmetry is 18.0 eV obtained by KLI which is compared 
with the HF (31.0 eV) and Slater (23.6) eV values. The 
overall conclusion can be drawn that the various orbital 
dependent potentials provide band structures similar to 
one another and all of them differ from HF. The largest 
deviation from experiment is found for the lowest valence 
band width (LVBW) which is rather narrow in the OEP 
methods. HF give too wide LVBW while Xa-LDA yields 
a value in closest agreement with experiment. 



those obtained by Hartree-Fock. The addition of the 
step-potential as a response part of the exact exchange 
(Eq. 13) to the Slater's potential considerably improves 
the quality of various physical properties calculated for 
polyethylene. To test the quality and reliability of our 
code we performed calculations in the atomic and molec- 
ular limit. The results can be brought into harmony with 
those obtained by others. 

In general, we find that the band structure calculated 
by different v x are similar with little qualitative differ- 
ence, while all of them differ significantly from HF which 
have valence bands much deeper and much higher virtual 
levels. The single-particle band gap obtained by the KLI 
or Slater approaches is too narrow when compared with 
the experiment. The shifted values are observed for the 
Fermi level, which differ from Hartree-Fock values signifi- 
cantly. The one electron energy levels are carefully exam- 
ined and an upward shift of eigenvalues observed when 
compared with those of HF, and also the counter shift of 
virtual levels. The deviation from the HF eigenvalues is 
more significant then in molecules or in atoms. In order 
to make a valuable comparison with Hartree-Fock, LDA 
and with the experimental gap, correction must be taken 
into account to the Kohn-Sham eigenvalue gap. The ex- 
change contribution to the derivative discontinuity is cal- 
culated by the KLI exchange potential. In this way the 
KS single particle gap is corrected by a single shot calcu- 
lation on top of exchange-only OEP. The exchange contri- 
bution to the band gap is somewhat less than twice the 
KS eigenvalue gap. According to the expectations the 
corrected gap is much closer to the Hartree-Fock value 
than the KS eigenvalue gap. 



Conclusions 

Density functional of the third generation treat both 
the kinetic and exchange energy exactly at the orbital de- 
pendent level of the theory. These orbitals come from a 
local potential and, due to the Hohenberg-Kohn theorem, 
are therefore functionals of the density. The optimized 
effective potential method, in the form given by Talman 
and Shadwick [1] allows one to treat orbital-dependent 
functionals in the framework of Kohn-Sham density func- 
tional theory. The OEP integral equations are difficult to 
solve for extended systems and therefore several approx- 
imations are introduced. Among them the one obtained 
by Krieger et al. (KLI) and another one derived by Grit- 
senko et al. have become known in the last few years. 

In this work we have applied the exchange-only KLI 
approximation to the OEP approach for insulating poly- 
mers. Calculations have been carried out for polyethy- 
lene with a new polymer code using the efficient KLI 
approximation to OEP. The results obtained by the 
"naked" Slater's potential and KLI are compared with 
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TABLE I. Comparison of our HF and KLI results 
(Clementi's double zeta basis) using the code DIOGENES 
[33] with values from the literature. 

Etot and e m are the total energy and the highest occupied 
orbital energy in a.u. 









KLI 01 "' 


KU other 


Be Etot 


-14.569 


-14.5730 a 


-14.5684 


-14.5723 b 




-0.3089 


-0.3093 


-.3092 


-.3089 


LiH Etot 


-7.954 


-7.987 b 


-7.954 


-7.987 c 




-7.8629 d 


-7.8620 e 








-0.2914 


-0.3017 b 


-0.3227 


-0.3011 c 




-.2870 d 


-.2857 e 







N 2 


Etot 


-108.892 ' 
-0.6352 f 


-108.994 d 
-0.6152 c 


-108.883 ' 
-0.6659 f 


-108.986 c 
-0.6818 c 


CH 2 


E tot 

Cm 


-38.8650 ' 
-0.3750 f 




-38.8632 ' 
-0.3595 f 







a [7] 
b [8] 
C [12] 



d with STO-3G basis 

e with STO-3G basis using Gaussian94 [43] program package 
f The large 11S/7P basis set is used for N2 and for CHq,. 



Fig. 1 The calculated valence band structure obtained by various 
exchange-only methods as a function of the dimensionless variable k 
with k = being the zone center and k = 1 the zone boundary. The 
Brillouin zone is that corresponding to a CH-z unit cell. Solid lines 
correspond to the occupied levels while the dashed lines to the virtual 
levels. The two lowest dashed curves of virtual levels are of particular 
interest (see text). The Clementi's double-^ basis set is used and 8 
neighbours are considered for the CH2 unit cell. 
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TABLE III. Calculated properties of polyethylene by vari- 
ous DFT methods using the Clementi's double-zeta basis set 
(9S/5P). 

Results with basis (11S/7P) are also given in the bottom of 
the Table. The sum of the calculated KS eigenvalue gap 
A^onint and the exchange derivative discontinuity A x is also 
given together with the exchange-only LDA and experimental 
gap (eV). The notations are the same as in Table I. 



TABLE II. Calculated properties of polyethylene by vari- 
ous DFT methods using the minimal basis set 
HF, Slater, SSP and KLI denote the Hartree-Fock, ex- 
change-only method with Slater's potential, Slater's poten- 
tial with orbital dependent step potential and the KLI ex- 
change-only OEP method. All the properties are in a.u. 
except the HOMO-LUMO gap which is given in eV. The 
Clementi's minimal basis set (7S/3P) is used. Neig gives the 
number of neighbouring cells taken into account. E tot , E x are 
the calculated total and exchange energies per CH2 unit cell. 
A x is the exchange component of the derivative discontinuity. 
The calculations are carried out for the CH 2 unit cell. 



neig 


HF 


Slater a 


SSP a 


KLI 


1 E tot (a.u.) 


-39.03535 


-39.01888 


-39.02768 


-39.02042 


E x (a.u.) 


-5.82128 


-5.84774 


-5.79080 


-5.7713 


e m (a.u.) 


-.6672 


-.5634 


-.4924 


-.3524 


e m+ i (a.u.) 


.3978 


.0074 


-.1084 


.2228 


gap (eV) 


29.0 


15.1 


16.4 


15.7 


A, (eV) 




-5.7 


-4.4 


-3.2 


5 


-38.88012 


-38.86891 


-38.87278 


-38.86981 




-5.75048 


-5.76397 


-5.7379 


-5.7297 




-.4747 


-.3524 


-.2312 


-.0871 




.3444 


-.0433 


.1310 


.2786 




22.3 


8.4 


9.9 


10.0 






-3.5 


-1.9 


-1.2 


8 


-38.88008 


-38.86890 


-38.87275 


-38.86954 




-5.75049 


-5.76400 


-5.7378 


-5.7284 




-.4760 


-.3518 


-.2318 


-.0972 




.3436 


-.0425 


.1303 


.275 




22.3 


8.4 


9.9 


10.1 






-3.5 


-1.9 


-1.1 


KS gap+Aj, 




4.9 


8.0 


9.0 



" [15] 



neig 


HF 


Slater a 


SSP a 


KLI 


5 Etot 


-39.01068 


-38.96683 


-38.99704 


-38.99311 


E x 


-5.8783 


-5.9368 


-5.8841 


-5.8405 




-.4049 


-.2951 


-.1628 


-.1023 


£m + l 


.1489 


-.1590 


-.0007 


.0284 


KS gap 


15.1 


3.7 


4.4 


3.6 


Ax 




3.4 


3.8 


1.9 


8 (9S/5P) 


-39.01056 


-38.96684 


-38.99760 


-38.99407 




-5.8780 


-5.9367 


-5.9036 


-5.8621 




-.3984 


-.2774 


-.1130 


-.0997 




.1558 


-.1655 


-.0755 


.0322 




15.1 


3.7 


4.4 


3.6 






3.2 


3.7 


2.4 


8 (11S/7P) 


-39.02326 


-38.99926 


-39.01361 


-39.00701 




-5.8977 


-5.9486 


-5.9036 


-5.8553 




-.3786 


-.2401 


-.1130 


-.0933 




.2320 


-.1151 


.0755 


.0958 




16.6 


3.4 


5.1 


5.1 






6.3 


7.6 


7.7 


KS gap+Az 




9.7 


12.7 


12.8 


LDA 








7.6 a 


EXP 








8.8 b 



a [15] 
b [24,38] 



TABLE IV. One-electron energies (a.u.) obtained by 
HF and KLI exchange-only approaches using the Clementi's 
(11S/7P) basis set and 8 neighbours for polyethylene (k = 1) 





HF moi 


KL jmo( 


Hppoty 


KLF oi » 


£1 


-11.289 


-10.201 


-11.219 


-9.820 


£2 


-.904 


-.774 


-.800 


-.511 


£3 


-.587 


-.537 


-.681 


-.401 


£4 


-.375 


-.359 


-.379 


-.093 


£5 


.056 


-.318 


.232 


.096 


£6 


.240 


-.058 


.232 


.109 


£7 


.314 


-.005 


.392 


.301 


£8 


.518 


.192 


.719 


.527 
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